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Theoretical predictions for the dynamic structure factor of a harmonically trapped Fermi su- 
perfluid near the BEC-BCS crossover are compared with recent Bragg spectroscopy measurements 
at large transferred momenta. The calculations are based on a random-phase (or time-dependent 
Hartree-Fock-Gorkov) approximation generalized to the strongly interacting regime. Excellent agree- 
ment with experimental spectra at low temperatures is obtained, with no free parameters. Theo- 
retical predictions for zero-temperature static structure factor are also found to agree well with the 
experimental results and independent theoretical calculations based on the exact Tan relations. The 
temperature dependence of the structure factors at unitarity is predicted. 

PACS numbers: 05.30.Fk, 03.75.Hh, 03.75.Ss, 67.85.-d 



Ultracold Fermi gases of 6 Li and 40 K atoms near Fes- 
hbach resonances provide a new paradigm for studying 
strongly correlated many-body systems At low tem- 
peratures, they display the intriguing crossover from a 
Bose-Einstein condensate (BEC) to a Bardeen-Cooper- 
Schrieffer (BCS) superfluid In the unitarity regime 
at the cusp of the crossover, a superfluid with nei- 
ther dominant bosonic nor fermionic character emerges 
that exhibits universal properties that might be found 
in other strongly interacting superfluids 0, 13], such as 
high-temperature superconductors or nuclear matter in 
neutron stars. This new superfluid has already been in- 
vestigated intensively [l], leading to several milestone 
observations, some of which still defy theoretical under- 
standing. Here we present a quantitative description of 
the recent two-photon Bragg spectroscopy measurement 
for this new superfluid [f|. 

Theoretical challenges in describing the BEC-BCS 
crossover arise from its strongly correlated nature: there 
is no small interaction parameter to set the accuracy of 
theories @. Significant progress has been made in de- 
veloping better quantum Monte Carlo simulations @- 
[To| and strong-coupling theories [il [n-14|, leading to 
the quantitative establishment of a number of proper- 
ties. These include equation of state fj, M. |l2l 11 5 — 
[l7| . frequency of collective oscillations la, ll°|. pair- 
ing gap (lol [20| , and superfluid transition temperature 
SEH- However, other fundamental properties, such as 
the single-particle spectral function measured by rf spec- 
troscopy [22,[2^| and the dynamic structure factor probed 
by Bragg spectroscopy Q, are not as well understood. 

In this Rapid Communication, we show that a random- 
phase approximation (RPA), generalized to the strongly 
interacting regime, is able to describe quantitatively 
the observed Bragg spectra for harmonically trapped 
6 Li atoms at large transferred momenta. This sur- 
prising result indicates that the RPA captures the es- 
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Figure 1: (color online). Quantitative comparison of the- 
oretical and experimental Bragg spectra [see Eq. ©]. The 
RPA prediction (lines) agrees well with the experimental data 
(empty squares) [|[ at the BEC-BCS crossover, with no free 
parameters. The spectrum is normalized so that the area be- 
low the curve is unity. The frequency is measured in units of 
the recoil energy of the atoms (see text). 



sential physics and constitutes a reasonable approxima- 
tion for the strongly interacting region of the BEC-BCS 
crossover, particularly the low temperature range ac- 
cessed by most experiments. The RPA method has pre- 
viously been used to study the dynamic structure factor 
(2^ and collective oscillations [25| of weakly interacting 
Fermi superfluids. A dynamic mean-field approach, iden- 
tical to the RPA but based on kinetic equations, was de- 
veloped to investigate structure factors [26| and collective 
modes [StJ of a uniform, strongly interacting Fermi gas. 
At finite temperatures, structure factors at the crossover 
were also studied using a pseudogap theory [2aj . 

Our main result is summarized in Fig. 1, which shows 
the normalized experimental Bragg spectra Q along with 
the RPA predictions. Excellent agreement is found, with 
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no free parameters. 

We begin by outlining briefly the RPA using the Hamil- 
tonian (hereafter h = 1), 



and the static structure factor is given by S a 



+U f dr^+(r)^+(r)^(r)^ t (r), 



Mr) 



(1) 



which describes a balanced spin- 1/2 {a =t,|) Fermi 
gas with mass M in a harmonic trap Vr(r), where 
fermions with unlike spins interact via a contact potential 
Uq8{y — r'). The total number of atoms N is tuned by the 
chemical potential fi and the bare interaction strength 
Uq is renormalized by the s-wave scattering length a, 
1/Uq + J2k A//k 2 = M/(47ro). In the superfluid phase, 
we treat the system as a gas of long-lived Bogoliubov 
quasiparticles interacting through a mean-field and con- 
sider its response to a weak external field of the form of 
5Ve i{<i-T-wt)_ The essential idea of the RPA is that there 

is a self-generated mean-field potential experienced by 
quasiparticles [2j|, associated with the local changes in 
the density distribution of the two spin species, SU = 
U J dr( (^V^ct Vv + Smip^ip^ + 5m*tp^), where 
5n a = 5n a (r, t) and 5m = 5m (r, t) are the normal and 
anomalous density fluctuations, respectively, which must 
be determined self-consistently. In the linear approxima- 
tion, the self-generated potential 5U plays the same role 
as the perturbation field when we calculate the dynamic 
response using a static BCS Hamiltonian as the reference 
system 24, 2a, 2f|. This leads to coupled equations for 
density fluctuations. The linear response is character- 
ized by a matrix consisting of all two-particle response 
functions: 

((«t«t)) ((%%)) {{n\rh)) ((n t m+)) 

((Hnt)) ((%™» ((™i™ + )> 

((mfif)) ((rhh^)) {(mm)) ((mm + )) 

((m + h^)) ((m + fii)) ((m + m)) ((m + m + )) 

where ((AB)) is the Fourier transform of the retarded 
function -iG (t - t') ([i(r, t), B(r', t')\). For simplicity, 
we abbreviate Xaa> = ((nafia')), X<ym = ((nam)), x<xm = 
((n CT m + )), Xmm = ((mm + )), and so on. By solving the 
coupled equations for density fluctuations, the standard 
RPA response function x can be ex pre ssed in terms of 
the static BCS response function x° |25fl . 



x = x° [i-Uox°G]~ 



(2) 



where Q — 5(r — i"')! ^ ® <r x ] is a direct product of two 
Pauli matrices Co and &x and the unit matrix 1 = i5(r — 
r')[<7o <8> <To]. The dynamic structure factor S acr i(oj) is 
related to the normal density response function by the 
fluctuation-dissipation theorem, 

S aa >(u) = T r — Imxaa' (u) , (3) 

7T [1 — exp {—lo/kbI )\ 



(2/iV) f^°° du>S aa ' (w). In the weak-coupling regime, 
Eq. @ can be solved by calcula ting x° f° r a thermal 
average of BCS quasiparticles [2~il . l25j 

Here, we extend the RPA to the strongly interacting 
regime with an arbitrarily large scattering length a, by 
properly renormalizing the bare interaction strength Uq 
and the two response functions Xmm an d X%,m-> wn ich 
was found to be suitable at the BEC-BCS crossover 

0, 

1 3 Oil . The ultraviolet divergence of these two functions 
[25| is canceled exactly by the small value of Uq, when 
the momentum cut-off goes to infinity. In homogeneous 
systems, a careful account of the divergent terms in the 
inverted matrix of the RPA equation @ leads to a concise 
expression for the response functions: 
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and 



xn = xtt - x?t + Xn, 



(x° n ) 2 , (4) 
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where the response functions with a tilde, i.e 
xlm + Ek M/k 2 — M/ (4™), become free from any ultra- 
violet divergence. The above equations were previously 
obtained by Combescot and collabrators using kinetic 
equations (see Eq. (B22) in Ref. Note that, we 

use a Leggett-BCS ground state without inclusion of the 
Hartree-Fock term in the quasiparticle spectrum. There- 
fore, in the BCS regime our treatment does not account 
for the leading interaction effect as in Rcfs. 0, 25 1. 



At the crossover, however, it does capture the dominant 
pairing gap. Note also that, the RPA method accounts 
for single particle-hole excitations. Higher correlations 
such as multi-particle-hole excitations are neglected. 

In the presence of a harmonic trap, the renormaliza- 
tion procedure becomes cumbersome because of the dis- 
crete energy levels. It is convenient to use a local den- 
sity approximation (LDA) that treats the system as a 
collection of many homogeneous cells with local chem- 
ical potential ^(r) = /i — Vr(r), where Vr(r) = 
M(u>^x 2 + uiyU 2 + w 2 z 2 )/2 is the harmonic trapping po- 
tential. The LDA treatment is valid for a large number 
of atoms such as N ~ 10 5 as in experiments. It has 
been used extensively in studying the static density pro- 
file of either atomic Fermi, Bose gases [3l| or Bose-Fermi 
mixtures [jj^. In the nuclear context, it has also been 
used to calculate the dynamic response function [33| . At 
a given temperature and scattering length, we solve the 
Leggett-BCS equation with local chemical potential for 
the local pairing gap and calculate the static response 
function x°, then solve the local RPA density response 
functions using Eqs. ^ and ([S]), and finally obtain the 
total RPA responses by integrating over the whole trap. 
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In our calculations, the interaction strength is charac- 
terized by the dimensionless parameter, l/(kpa), where 
k F = ^2ME F is the Fermi wave vector and the Fermi 
energy is Ep = {'iNu! x uiyijj z ) 1 /' i . 




Figure 2: (color online). Zero temperature spin parallel 
Sff(q,Lu) (dashed lines), anti- parallel S^_i(q,to) (dot-dashed 
lines), and total dynamic structure factor S^q, to) = 2[S^ + 
S^i] (solid lines) across the BEC-BCS crossover: l/k F a = 0.5 
(a), 0.0 (b), and —0.5 (c). The negative weight in Sfi at 
about the recoil energy is consistent with the exact sum rule 
/ u!Sfi(q, u>)du! — [281 ]. 

Figure 2 shows the zero-temperature spin parallel, 
anti-parallel, and total dynamic structure factor at a 
transferred wave-vector q = hk F in the BEC-BCS 
crossover, calculated using the above RPA procedure for 
a trapped Fermi gas. In addition to a broad response 
at the recoil energy lur = q 2 /(2M) = 25E F caused by 
resonant scattering of atoms, a much narrower peak de- 
velops at about lor/2 with increased coupling. The peak, 
commonly referred to as the quasielastic peak in the liter- 
ature, is found by the recent theoretical calculation [26[ 
and the observation of Bragg spectroscopy This is 
simply the Bogoliubov- Anderson phonon mode of a Fermi 
superfluid at large wave-vectors, which evolves continu- 
ously into a Bogoliubov mode of molecules towards the 
BEC limit (2rj j. The molecular peak is mostly evident 
in Sf± as there is no background atomic response. Mea- 
surments of ma y a lso help establish the presence of 
Fermi superfluidity [28|. 

To make a quantitative comparison with the experi- 
mental spectra, we calculate the momentum imparted to 
the Fermi cloud, the quantity measured directly in the 
Bragg scattering experiment |5l. |34|: 



1 r°° 

P(q,w) cx — / dw'S(q,u/' 

7TCT J_ ao 



smc 



(6) 



where sinc(x) = sm(x)/x and the energy resolution 
a = 2 Jt Br is set by the experimental Bragg pulse du- 
ration (T Br = 40 ^s) H. We find a w 0.68E F 
0.027wij. Figure 1 presents a comparison of the exper- 
imental data (open squares) with the RPA predictions 
(lines) for the Bragg spectra normalized in such a way 
that / V(q,uj)duj = 1. With no free parameters, our 



RPA predictions agree well with the experimental results 
in the unitarity regime (1/kpa = 0.0 and 0.2) and BEC 
regime {1/kpa = 0.8). The agreement on the BCS side 
(1/kpa = —0.5), however, becomes worse. The quantita- 
tive agreement around unitarity is very compelling, since 
the RPA was assumed to be unreliable in the (strongly 
interacting) regime of large pair fluctuations. Our com- 
parison indicates that the RPA is able to describe the dy- 
namical properties of the BEC-BCS crossover, at least at 
zero temperature and large momenta. High order multi- 
particle-hole excitations, absent in the RPA theory, seems 
to be negligibly small at large momenta. More studies 
are needed to understand this. Note finally that, the 
somewhat poorer agreement at 1/kpa = —0.5 can be at- 
tributed to the mean-field shift, which is ignored in the 
RPA but dominates for sufficiently weak interactions. 
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Figure 3: (color online). Quantitative comparison between 
theory and experiment for the zero temperature static struc- 
ture factor across the crossover. For Sfi, with no free pa- 
rameters our RPA prediction (plus symbols) agrees well with 
the experimental data for 5(q, uj) — 1 (solid circles with error 
bars) [35fl and an independent theoretical result based on the 
exact Tan relations (solid line) (36| . At large transferred mo- 
mentum, S 1 ^ ~ 1. The inset highlights the RPA prediction 
with respect to the Tan-relation result in the BCS regime. 

The agreement between the RPA theory and the Bragg 
experiment is further confirmed by comparing the spin 
anti-parallel static structure factor at zero temperature, 
as reported in Fig. 3. Experimentally, the static struc- 
ture factor can be measured model-independently by in- 
voking the /-sum rule (3||; while, theoretically, it can be 
determined very accurately using the exact Tan relations 
and the known equation of state [36[. It is evident from 
Fig. 3 that the RPA prediction fits very well with the 
experimental data, as well as the independent theoreti- 
cal result based on the Tan relations. In particular, the 
two theoretical predictions are nearly indistinguishable 
on the BEC side with 1/kpa > 0. However, they differ 
towards the BCS limit, as highlighted in the inset. The 
discrepancy is consistent with Fig. Id where the RPA 
predicts less pairing and hence lower S(q). 
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Figure 4: (color online). Temperature dependence of the 
dynamic (a) and static (b) structure factor for a unitary 
Fermi gas in harmonic traps at q = 5fci?. According to 
the Tan relation, 5 n (q) ^ 128C/[175£ 1/4 (q/fcF)] at T = 
[35| . where £ and £ are the universal parameters at unitar- 
ity [l7| . The symbols in (b) show predictions using theo- 
retically or experimentally determined £ and f: ENS exper- 
iment (square) [l7| . Gaussian pair fluctuation theory (cir- 
cle) fl2fl . and self-consistnent theory (triangle) [37l |. Here, 
T c ~ 0.37T F = 0.37E F /k B . 



A more stringent test of the RPA theory may be pro- 
vided by the temperature or momentum dependence of 
dynamic and static structure factors. In Fig. 4, we pre- 
dict the dynamic and static structure factor as a func- 
tion of temperature for a trapped Fermi gas at unitar- 
ity, which will be investigated in future experiments. As 
anticipated, the pair (atomic) response increases (de- 
creases) with decreasing the temperature, leading to a 
monotonic decay of the static structure factor. 

The present RPA theory is most likely valid only in a 
narrow temperature window near T = 0. With increas- 
ing temperature, the pairing gap decreases and thermal 
pair fluctuations increase. The RPA will eventually break 
down at a characteristic temperature Trpa(% T c ). This 
is evident in Fig. 4b where the spin anti-parallel static 
structure factor vanishes unphysically above the super- 
fluid transition temperature. 

At low transferred momenta, quantum fluctuations are 
likely to increase and the RPA theory will become less 
reliable. To overcome these limitations, we could use 
a Cooperon-mediated interaction (many-body T-matrix) 
to replace the bare contact interaction [38[, or use the 
phenomenological Landau parameters for the mean-field 
shift |39|| . as determined from thermodynamic measure- 
ments |I6| or quantum Monte Carlo simulations. 

In summary, we have used a strong-coupling RPA the- 
ory to calculate the dynamic and static structure factors 
of a trapped Fermi gas at the BEC-BCS crossover. The 
theory is quantitatively applicable at low temperatures 
and large transferred momenta, as confirmed by the ex- 
cellent agreement with the experimental Bragg spectra. 
The RPA theory thus seems to provide a novel starting 
point for investigating dynamic properties of a strongly 



interacting Fermi gas at finite temperatures and low mo- 
menta. 
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